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We search for CP violation in Cabibbo-suppressed charged D meson decays by measuring the dif- 
ference between the CP violating asymmetries for the Cabibbo-suppressed decays — ^ K'^K'n'^ 
and the Cabibbo-favored decays Df — > K^K~n^ in the K^K~ mass region of the resonance. Us- 
ing 955 fb~^ of data collected with the Belle detector we obtain Agp^*^"* = (-|-0.51±0.28±0.05)%. 
The measurement improves the sensitivity of previous searches by more than a factor of five. We 
find no evidence for direct CP violation. 

PACS numbers: 13.25.Ft, ll.30.Er, 14.40.Lb 



Studying CP asymmetries in D meson decays provides 
a promising opportunity to search for new physics (NP) 
beyond the Standard Model (SM) [1]. Here we study 
CP asymmetries in charged Z3+ — > (pn^ and Df — > (pn^ 
decays 01. The observable of interest is 
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where F is the partial decay width. This time-integrated 
asymmetry arises from CP violation (CPV) in decay am- 
plitudes. Within the SM, CPV in D decay amplitudes is 
predicted to be very small. The largest effect occurs for 
singly Cabibbo-suppressed (SCS) decays such as £>+ 

, which are governed by the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements VcsV*,,. However, even 
for these decays Acp is predicted to be only 0(0.1%) or 
less [Sj. In contrast, several NP models predict Acp to 
be as large as 1%. Experimentally, to cancel detector- 
induced asymmetries and other systematic effects, we 
measure the difference 



N{D+)- N{D-) _ N{D+) - N(Dj) 
N{D+) + N{D-) ' NiDt) + NiDs) ' 
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where the second term corresponds to the Cabibbo- 
favored (CF) decay Df — 7> 07r+. This CF decay is gov- 
erned by the CKM matrix elements V^^V*^ and is ex- 
pected to have negligible Acp Q; thus, a measurement 



of A^rcc probes A^p~^'^^ . Measuring a relatively large 



value would be interpreted as evidence for NP. Previ- 
ously, CPV in D meson SCS decays has been searched 
for in several final states @. No significant asymme- 
tries were found, with the best sensitivities rangingfrom 
0(0.2%) to 0(2%) depending on the decay mode g-Q. 

The measurement is based on 955 fb"-'^ of data recorded 
with the Belle detector [§| at the KEKB asymmetric- 
energy e'^e~ collider [lo| . which primarily operated at 
the center-of-mass (CM) energy of the T(4S') resonance 
and 60 MeV below. A fraction of the data was recorded at 
the T(IS'), T(25), T{3S), and T(5S') resonances; these 
data are included in the measurement. The Belle detec- 
tor is described in detail elsewhere 0, [llj: in particu- 
lar, it includes a silicon vertex detector (SVD), a central 
drift chamber, an array of aerogel Cherenkov counters, 
time-of-flight scintillation counters, an electromagnetic 
calorimeter and a muon detector. 

We reconstruct the decays Z?^^ — 07r+ in the (p — )■ 
K^K^ decay mode. Each final state charged particle 
is required to have at least two associated SVD hits in 
each of the two measured coordinates. To select pion and 
kaon candidates, we impose standard particle identifica- 
tion criteria [13]. The identification efficiencies and the 
misidentification probabilities are about 90% and 5%, re- 
spectively. In addition, we require loose proton veto cri- 
teria for kaon candidates and loose lepton veto criteria 
for pion candidates, since we found that a considerable 
fraction of background (23%) involves misidentified pro- 
tons and leptons. D meson daughter particles are refitted 
to a common vertex and the D meson candidate is con- 
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strained to originate from the e+e" interaction region. 
Confidence levels exceeding fO~'^ are required for both 
fits. In order to reject D mesons produced in B meson 
decays, the D meson momentum in the e'^e~ CM system 
must satisfy > 2.5 GeV/c, for the data taken below 
T(5S'), and p*^ > 3.1 GeV/c for the T(5S) data. 

We accept candidates in the invariant mass regions of 
D and (I) mesons, 1.80 GeV/c^ < Mkk^ < 2.05 GeV/c^ 
and Mkk < 1-07 GeV/c^. For the small fraction of 
events with multiple candidates (4.6%), we select a sin- 
gle best candidate: the one with the smallest of the 
production and decay vertex fits. 

We study background using a generic Monte 
Carlo (MC}_simulation based on EVTGEN and 
GEANT3 We find that the main component (97%) 
is the combinatorial background whose shape in Mkktv 
can be fitted well with an exponential function. Other 
background components are mainly due to decays of 
charm particles and have a complicated structure in 
Mkk-k- However, their fractions are sufficiently small 
that the structure is obscured by the statistical fluctu- 
ations of the main background component. The combi- 
natorial background can be further divided into random 
combinations of a correctly reconstructed (j) meson and a 
7r+ (42%), and the rest (58%). 

To improve the purity of the and Z?+ data sample 
we require: \Mkk — rn^] < 16 MeV/c^, where is 
the nominal mass of 0, p^r > 0.38 GeV/c, where p-y^ is 
the laboratory momentum of 7r+, and | cos^heil > 0.28, 
where 9hei is the angle between K~ and momenta in 
the rest frame of (j). These selection criteria are obtained 
by minimizing the expected statistical error on AArec 
using signal and background samples from the generic 
MC simulation. The simulation has been tuned prior 
to running the optimization procedure to reproduce the 
mass resolutions of the D signals in data and the signal- 
to-background ratios of the data. 

The measured asymmetry Aj-ec can be written as the 
sum of several contributions that are assumed to be 
small: 



Arcc ^ Acp + Afb {cose*) + A, 



KK 
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In addition to the intrinsic asymmetry, Acp^ there is a 
forward-backward asymmetry (Afb) in the production 
of D mesons in e~^e~ — cc arising from 7 — Z*^ interfer- 
ence and higher-order QED effects. This term is an odd 
function of the cosine of the D meson production po- 
lar angle 6* in the CM system and could differ between 
and due to fragmentation effects. Furthermore, 
there are contributions due to asymmetry in the recon- 
struction efficiencies of oppositely charged kaons (A^^) 
and pions (A^). The term ^ = for ^- K+K' de- 
cays. However, the interference with other intermediate 
states in the decay D^^-^ K^K~it^ introduces a small 
difference in momentum distributions of the same-sign 



same sign 
opposite sign 




FIG. 1: Background subtracted momentum distributions of 
kaons of the same and opposite charges relative to that of 
the D meson: (a) for decays and (b) for D^t decays. 
Background is taken from sidebands in Mkktv- 



and the opposite-sign kaons (where the sign is relative to 
the Dj-^-j meson charge), as shown in Fig. [TJ The differ- 
ence in momentum distributions in combination with the 
kaon detection asymmetry A^ leads to a non-zero A^^ , 
as explained below. 

We define the intrinsic laboratory phase space distribu- 
tion of the kaon pair by P(a;i, 3:2), where = (pi,cos0i) 
and X2 = (p2,cos02) label the phase space variables of 
same-sign and opposite-sign kaons, respectively. The 
measured same-sign and opposite-sign single kaon dis- 
tributions are obtained from the intrinsic P(a;i,a;2) by: 

Pi2{xi2)= '^^^^^'^^^ a;2)e(a:2(i))rfa;2(i) 

where e{x) is the phase-space dependent detection effi- 
ciency. The numbers of detected positively and nega- 
tively charged D mesons are 



j J dxidx2P{xi,X2)eK±{xi)eKT{x2), 



(5) 



where e^±(a:) = e(a;)(l ± A^{x)) are the efficiencies of 
the as functions of kaon phase space x. From this 
and neglecting terms quadratic in A^ one obtains: 



A^''= l{Pi{x)-P2{x))Af{x)dx, 



(6) 



where Pi{x) and Piix) are normalized distributions of 
the detected same-sign and opposite-sign kaons, respec- 
tively, given by Eq. ^ and the integration runs over the 
kaon phase space x = {p,cos6). 

The last term in Eq. ([3]) is a function of pion mo- 
mentum and polar angle in the laboratory frame. In 
the difference of measured and Z?+ asymmetries, 
provided the measurement is done in bins of the three- 
dimensional (3D) phase space (cos^*,p7r,cos07r), the last 
term in Eq. ^ cancels. 
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In the above equation we assume that the intrinsic 
A^p'^'^^ is neghgible, as discussed in the introduction. 
We use 10x10x10 equal size bins of the 3D phase space 
with < 5 GeV/c. The yields of D+ , D- , D+ and 
Dj decays are determined from a binned likelihood fit 
to the MxK-n distributions in each sufficiently populated 
3D bin. We require at least 100 entries in the histogram 
in order to perform the fit. To parameterize the non- 
Gaussian signal shape with as few parameters as possi- 
ble, we use the distribution of pulls determined with MC 
simulation; the pulls are calculated as {Mkktt ~T^)/o'rm 
where m and a^, are the 3D bin dependent mean and 
standard deviation of the or Df invariant mass dis- 
tributions. The pull distribution is fitted with a sum of 
four Gaussians to obtain their fractions /f"", mean po- 
sitions x^"" and the widths erf"". The signal shape for 
the decays with no final state radiation (FSR) is param- 
eterized with: 



S4g{x) 
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where st = o-f^'V and = a;f""(T?"" -f xq. The normal- 
ized shape given by Eq. ([8|) has two free varying param- 
eters: the position xq and the width a. 

The pull distribution is found to be cos 9* dependent; it 
becomes asymmetric at the edges of cos 9* space. Thus, 
to improve the mass fits, we use cosO* dependent pull 
parameters /f"", a;f"" and erf"" in Eq. (|8]), which are 
obtained from fits to simulated distributions in ten bins 
of cos 9* . 

The^ape of the FSR tail 5'FSR(a;) and its fraction 
Pfsr [15i are taken from MC simulation using PHO- 
TOS [16| to simulate FSR. The normalized signal shape 
is then S{x) = (1 - PFSRjSigix) + PFSR5'FSR(a;)- The 
background is parameterized with an exponential func- 
tion of a quadratic polynomial, B{x) = g"-+bx+cx ^ 
free parameters a, b, and c. 

The parameterization of the Mkktt distribution in- 
cludes two signal peaks and the background, and has 
eight free varying parameters: two yields {No+, ^dJ, 
two peak positions {m£i+, m^)^), the width of £)+ peak 
(cr£)+) and three background parameters (a, b, c). The 
ratio of the and peak widths, / = crD,/o'D+7 is 
fixed from MC simulation in order to ensure stable fit- 
ting. 

Of the approximately 700 sufficiently populated invari- 
ant mass distributions (-D^j and -Dj~j per each 3D bin), 
658 are fitted successfully. The quality of fits is good: 
the mean of the normalized distribution is 1.000 and 
the r.m.s is 0.090 for 242 degrees of freedom; the corre- 
sponding confidence level distribution is uniform. From 
the fitted yields in 3D bins we calculate the and Df 
asymmetries, and the asymmetry differences AA^-cc- We 
consider only those bins in which the yield has a signifi- 
cance greater than 3(t; this requirement must be fulfilled 
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FIG. 2: Sum of invariant mass distributions with the sum of 
fitted functions superimposed for (a) positively and (b) nega- 
tively charged D mesons. The plots beneath the distributions 
show the residuals. 



for all four measured yields in a bin. The sum of mass 
distributions in these bins is shown in Fig. [5J We find 
237525 ± 577 and 722871 ± 931 decays. The 
residuals of the sum of all the successfully fitted distri- 
butions do not show any significant structure. 

The asymmetry difference in each bin is then corrected 
with AAf ^ for that bin: 



rcc -^rcc 



AA 



KK 



(9) 



The corrections are determined using Eq. © and the ex- 
perimental data for the -Pi (a;) and P^^^) distributions. In 
particular, we calculate A^^ for events in the signal win- 
dow, rriB ± 15 MeV/c^, and subtract the asymmetry for 
events in an equal width sideband displaced ±20 MeV/c^ 
from the nominal D meson mass tod. 

The kaon asymmetry A^ ^ which is needed in Eq. ([S]), 
is measured using — > 07r+ and — > X^tt"*" de- 
cays; for the latter decay the measured asymmetry can 
be expressed as Aj-oc = Acp + ApB — A^ + A'^ . By as- 
suming negligible CP violation (both are CF decays) and 
the same forward-backward asymmetry, and by neglect- 
ing the A^^ term in Eq. ([3]), the difference of measured 
asymmetries is equal to A^ 17]. 

The procedure is similar to that used in Q. First we 
determine the asymmetry of Df — > (f)Tr~^ in 3D bins using 
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FIG. 3: CP-violating asymmetry (a) and forward-backward 
asymmetry difference (b) in bins of |cosS*|. The horizontal 
line in (a) is a constant fit to the data points. 



the fitted yields. This asymmetry map is used to weight 
— ^ iir~7r+ events in order to determine the /D 
corrected yields in bins of the kaon phase space. The 
yields are obtained by the sideband subtraction method 
as in The kaon asymmetry map is then calculated 

from the corrected /D yields in the range < pk < 
4 GeV/c and — 1 < cos 9k < 1 divided into 10x10 equal 
size bins. 

The following corrections are obtained with Eq. (jG)) for 
the total 3D phase space: ^ = (-h0.060± 0.013)% for 
D+, Af^ = (-0.051 ± 0.012)% for D+ and AAf ^ = 
(-1-0.111 ± 0.025)%; the difference is not zero because of 
the opposite signs of the momentum asymmetries in 
and decays, as shown in Fig.[T] The uncertainties are 
due to statistical variations of and Pi {x) — {x) ; the 
error on ^Af^ is included in the systematic uncertainty. 

The corrected asymmetry differences in 3D bins, A^J;™ 
defined by Eq. are used to calculate error- weighted 
averages in bins of cos 6* ; error- weighted averages are ob- 
tained with the least squared fit. Finally, A^p'^'^^ and 
/S.Afb are extracted by adding/subtracting the asymme- 
try difference in opposite bins of cos 9* 



^A'i°l{cos9*) + ^A'iZ{-cos6*) ^ ^^^^ 
^^^^^ _ AAgg,-(cosr)-AAi:g^(-cosr) ^^^^ 
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The results are shown in Fig.[3l By fitting the data points 
of Fig. [3^ with a constant we obtain A^p^"^"^ = (0.51 ± 
0.28)%, where the error is statistical only. The result is 
consistent with zero within 1.8 standard deviations. 

Figure Wp shows the difference in forward-backward 
asymmetries. The test with respect to /^ApB = 
gives x^/ndf = 10.57/5, which corresponds to a confi- 
dence level of 6%; no significant difference is found be- 
tween the forward-backward asymmetries for and 
A fit to a constant yields a value of /S.Afb — 
(0.25±0.28(stot.))%. 

We consider five significant sources of systematic un- 
certainties (Table HI). As discussed before, the Af^ cor- 



TABLE I: Summary of systematic uncertainties in A^^p 



Source 


Uncertainty (%) 


Af^ corrections 


0.025 


3D binning 


0.026 


Invariant mass binning 


0.022 


Fitting procedure 


0.018 


Selection of fit results 


0.020 


Sum in quadrature 


0.050 



rections are uncertain to 0.025%. The impact of 3D bin- 
ning is studied by changing the binning from 10x10x10 
bins to 20x10x10, 10x20x10 and 10x10x20 bins; we 
obtain a 0.026% variation in A^p^'^'^ . By doubling 
the number of bins in the invariant mass histograms a 
variation of 0.022% is obtained. The impact of signal pa- 
rameterization is studied by replacing the four-Gaussian 
shape with a triple Gaussian shape and the impact of 
background parameterization is studied by replacing the 
default parametrization with a simple exponential func- 
tion. We also vary the range in which we fit the distri- 
butions; all these changes give a 0.013% variation in the 
result. The uncertainty of the width ratio /, which is 
fixed in the fit, propagates into a 0.012% uncertainty in 
the result. By adding the last two numbers in quadrature 
we obtain an estimate of 0.018% for the systematic un- 
certainty of the fitting procedure. The last source is the 
selection of fit results; by changing the requirement from 
N/aN > 3 to N/aN > 5 we obtain a 0.020% variation in 
the result. We estimate the total systematic uncertainty 
by summing individual contributions in quadrature; we 
obtain 0.050%. 

In summary, we searched for CP violation in the de- 
cays D"*" — )■ 4nr~^ by measuring the CP violating asym- 
metry difference between Cabibbo-suppressed {D~^) and 
Cabibbo- favored (-D^) decays in a mass region around 
the <f> resonance, ± 16 MeV/c^. We have made no at- 
tempt to disentangle the (j) from other intermediate reso- 
nances in this mass region. Using 955 fb^^ of experimen- 
tal data collected with the Belle detector and assuming 
negligible CPV in CF decays we measure: 



A 



CP 



= (+0.51 ±0.28 ±0.05)%. 



(12) 



The result shows no evidence for CP violation and agrees 
with SM predictions. Previously, the most precise mea- 
surements were from CLEO [l^ and BaBar J^] ; our mea- 
surement improves the precision by more than a factor 
of five. We also measure for the first time the differ- 
ence in the forward-backward asymmetries of I?+ and 
Df mesons and find no significant deviation from zero. 
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